On inverse problems for the multidimensional 
relativistic Newton equation at fixed energy 



Alexandre Jollivet 



Abstract. In this paper, we consider inverse scattering and inverse boundary value 
problems at sufficiently large and fixed energy for the multidimensional relativistic 
Newton equation with an external potential V , V G C 2 . Using known results, we 
obtain, in particular, theorems of uniqueness. 
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1. Introduction 
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1.1 Relativistic Newton equation. Consider the Newton equation in the relativistic 
case (that is the Newton-Einstein equation) in an open subset of K n , n > 2, 
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where V G C 2 (0,IR) (i.e. there exists V G C 2 (M n ,lR) such that V restricted to 
^ is equal to V) and x = x(t) is a C 1 function with values in 0. 

By || V ||c 2 we denote the supremum of the set {|5^T^(x)| | x G 0, j = 

(ii,..,in)G(Nu{o}) w ,Er=iii<2} ; 

The equation (1.1) is the equation of motion of a relativistic particle of mass 
m = 1 in an external scalar potential V (see [E] and, for example, Section 17 of 
[LL]). The potential V can be, for example, an electric potential or a gravitational 
potential. In this equation x is the position of the particle, p is its impulse, t is the 
time and c is the speed of light. 

For the equation (1.1) the energy 



E = ^l + ^- + V(x(t)) 

is an integral of motion. We denote by B c the euclidean open ball whose radius is 
c and whose centre is 0. 

In this paper we consider the equation (1.1) in two situations. We study 
equation (1.1) when 



O = D where D is a bounded strictly convex open subset of M n , n > 2, 
with C 2 boundary. 
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And we study equation (1.1) when 

(1.2b) = R n and \d J x V(x)\ < (3^(1 + \x\)- a -^ , x G M. n , 

for \j\ < 2 and some a > 1 (here j is the multiindex j G (N U {0}) n , |j| = XliLi ii 
and are positive real constants). 

For the equation (1.1) under condition (1.2a), we consider boundary data. 
For equation (1.1) under condition (1.2b), we consider scattering data. 

1.2 Boundary data. For the equation (1.1) under condition (1.2a), one can prove 
that at sufficiently large energy E (i.e. E > E(\\V\\c 2 , D)), the solutions x of 
energy E have the following properties (see Subsections 3.1, 3.2 and 3.3 of Section 
3): 

for each solution x(t) there are t\,t2 G K, t\ < t2, such that 
(1.3) x G C 3 ([*i,t 2 ],R n ),a:(ti),x(t2) e G D for t e]t u t 2 [, 

x(si) ^ x(s 2 ) for si,s 2 G [ti,t 2 ],si 7^ s 2! 



(1 4) c ^ s ^ nc ^ P°^ n ^ s Go, Q £ D, there is one and only one solution 

x(t) = x(t, E, q , q) such that x(0) = qo, x(s) = q for some s > 0. 

Let (qo,q) be two distinct points of 3D. By s(E,qo,q) we denote the time 
at which x(t, E,q ,q) reaches q. By k(E,q ,q) we denote the velocity vector 
x(s(E, q , q), E, qo, q). We consider k(E, q , q), q , q G dD, q ^ q, as the boundary 
value data. 

1.3 Scattering data. For the equation (1.1) under condition (1.2b), the following 
is valid (see [Y]): for any (y-,x-) G B c x IR n , v_ 7^ 0, the equation (1.1) has a 
unique solution x G C 2 (R, lR n ) such that 

(1.5) x(t) = + +y-(t), 

where y~(t) — > 0, ?/-(£) — > 0, as t — > —00; in addition for almost any (i>_,x_) G 

(1.6) x(t) = v+t + x + +y+(t), 

where v + 7^ 0, \v+\ < c, v+ = a(t>_,x_), x+ = b(v-,x-), y+(t) — > 0, y+(t) — > 
0, as £ — > +00. 

For an energy £ > c 2 , the map 5 E : §e x R" ^ § E x R n (where § E = {v G 
-B c I M = cy^ 1 — (^) }) given by the formulas 

(1.7) v+ = a(v-, X-), x+ = b(v-,X-), 
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is called the scattering map at fixed energy E for the equation (1.1) under condition 
(1.2b). By D(Se) we denote the domain of definition of Se- The data a(v-,X-), 
b(v-,X-) for (v-,X-) G V(Se) are called the scattering data at fixed energy E 
for the equation (1.1) under condition (1.2b). 

1.4 Inverse scattering and boundary value problems. In the present paper, we con- 
sider the following inverse boundary value problem at fixed energy for the equation 
(1.1) under condition (1.2a): 

Problem 1 : given k(E, q , q) for all (qo, q) G 3D x dD, 

qo 7^ q, at fixed sufficiently large energy E, find V. 

The main results of the present work include the following theorem of uniqueness 
for Problem 1. 

Theorem 1.1. At fixed E > E(\\V\\c^, D), the boundary data k(E,q,q ), 
(qo,q) G dD x dD, qo ^ q, uniquely determine V . 

Theorem 1.1 follows from a reduction of Problem 1 to the problem of deter- 
mining an isotropic Riemannian metric from its hodograph and from Theorem 3.1 
(see Section 3). 

In the present paper, we also consider the following inverse scattering problem 
at fixed energy for the equation (1.1) under condition (1.2b): 

Problem 2 : given Se at fixed energy E, find V. 

The main results of the present work include the following theorem of uniqueness 
for Problem 2. 

Theorem 1.2. Let A G IR + and let D be a bounded strictly convex open 
subset ofR n withC 2 boundary. LetV 1 ,V 2 G Cg(M n ,M), max(|| Vi || C 2, H^Hc 2 ) < X 
and supp(Vi) U supp(V2) C D. Let S l E be the scattering map at fixed energy E 
subordinate to Vi for i = 1,2. There exists a nonnegative real constant E(X,D) 
such that for any E > E(X, D), V\ = V 2 if and only if S E = S E . 

Theorem 1.2 follows from Theorem 1.1 and Proposition 2.1. 

Remark 1.1. Note that for V G Cg(M n ,M), if E < c 2 + sup{V(x) \ x G M. n } 
then Se does not determine uniquely V. 

Note also that reducing Problem 1 to the problem of determining an isotropic 
Riemannian metric from its hodograph, one can give also stability estimates for 
Problem 1 under the assumptions of Theorem 1.1. 

1.5 Historical remarks. An inverse boundary value problem at fixed energy and at 
high energies was studied in [GN] for the multidimensional nonrelativistic Newton 
equation in a bounded open strictly convex domain. In [GN] results of uniqueness 
and stability for the inverse boundary value problem at fixed energy are derived 
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from results for the problem of determining an isotropic Riemannian metric from 
its hodograph (for this geometrical problem, see [MR], [B] and [BG]). 

Novikov [N2] studied inverse scattering for nonrelativistic multidimensional 
Newton equation. Novikov [N2] gave, in particular, a connection between the in- 
verse scattering problem at fixed energy and Gerver-Nadirashvili's inverse bound- 
ary value problem at fixed energy. Theorem 1.2 of the present work is a general- 
ization of theorem 5.2 of [N2] to the relativistic case. 

Inverse scattering at high energies for the relativistic multidimensional New- 
ton equation was studied by the author (see [Jl], [J2]). 

As regards analogs of Theorems 1.1, 1.2 and Proposition 2.1 for nonrelativistic 
quantum mechanics see [Nl] , [NSU] , [N3] and further references therein. As regards 
an analog of Theorem 1.2 for relativistic quantum mechanics see [I]. As regards 
results given in the literature on inverse scattering in quantum mechanics at high 
energy limit see references given in [J2]. 

1.6 Structure of the paper. The paper is organized as follows. In Section 2, we give 
some properties of boundary data and scattering data and we connect the inverse 
scattering problem at fixed energy to the inverse boundary value problem at fixed 
energy. In Section 3, we obtain a theorem of uniqueness and stability (Theorem 
3.1) for the inverse boundary value problem. This theorem is a generalization to 
relativistic case of theorem 4 of [GN]. 

Acknowledgement. This work was fulfilled in the framework of Ph. D. thesis 
researchs under the direction of R.G. Novikov. 



2. Scattering data and boundary value data. 

2.1 Properties of the boundary value data. Let D be a bounded strictly convex 
open subset of M n , n > 2, with C 2 boundary. 

At fixed sufficiently large E (i.e. E > E(\\V\\ C *,D) > c 2 + sup xeD V(x)) solu- 
tions x(t) of the equation (1.1) under condition (1.2a) have the following properties 
(see Subsections 3.1, 3.2 and 3.3 of Section 3): 

for each solution x(i) there are t±, t 2 G K, t\ < t 2 , such that 
x e C 3 ([t 1 ,t 2 ],M n ),x(t 1 ),x(t 2 ) e dD,x(t) e D for t e]t l ,t 2 [, 

(2.1) x(si) ^ x(s 2 ) for si, s 2 G [t u t 2 ], Sl ^ s 2 , x(*i)iV(x(ti)) < 
and x(t 2 )N(x(t 2 )) > 0, where N(x(ti)) is the unit outward 
normal vector of dD at x(U) for % = 1,2; 

for any two points qo,q G D,q ^ q , there is one and only one solution 

(2.2) x(t) = x{t, E, q , q) such that x(0) = qo,x(s) = q for some s > 0; 
x(0,E,q ,q) G C X {{D x D)\G,R n ), where G is the diagonal in D x D, 
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(where by u x(0, E, q , q) G C l ((D x D)\G, R n )" we mean that there exists an open 
neighborhood Q of D such that x(0, E, qo, q) is the restriction to (D x D)\G of a 
function which belongs to C 1 ((Q x 0)\A) where A is the diagonal of O x O). Let 
E > E(\\V\\c*, D). Consider the solution x(t,E,q ,q) from (2.2) for q ,q G dD, 
qo j£ q. We remind that s = s(E, qo, q) is the root of the equation 

x(s,E,q ,q) = q, s > 0, 

and we remind that k(E,qo,q) = x(s(E,qo,q),E,qo,q). We consider k(E,qo,q), 
qo, q G dD, qo 7^ q as the boundary value data. 
Let ko(E, qo, q) = x(0, E, q , q). Note that 

k (E,q ,q) = - k(E,q,q ), 



(2 ' 3) iwit M /1 ^-^0) 

|fco(^,?o,9)| =c\ 1 - 



c 2 



for £ > £(||y|| C 2,L>) and (g, g ) e (9£> x dD)\dG. 

2.2 Boundary data for the non relativistic case. If one considers the nonrelativis- 
tic Newton equation in D instead of the equation (1.1) under condition (1.2a), 
one obtains the existence of a constant £"(||y||c 2 , D) such that the solutions x(t) 
of the nonrelativistic Newton equation with energy E = ^\x(t)\ 2 + V(x(t)), E > 
E'(\\V\\c2, D), also have properties (2.1) and (2.2) (see [GN]). Hence one can define 
the time s'(E,q ,q) and the vector k'(E,q ,q) for E > E'(\\V\\c*, D), (qo,q) G 
(dD x dD)\dG, as were defined s(E,q ,q), k(E,q ,q) for E > E(\\V\\ C »,D), 
(q ,q) G (dD x dD)\dG. In [GN], s'(E,q ,q) and k'(E,q ,q) for E > E'(\\V\\ C 2, 
D), (qo,q) G (dD x dD)\dG, were taken as boundary value data for the multidi- 
mensional nonrelativistic Newton equation and [GN] obtains, in particular, that 
s'(E, q , q) given for all E > E'(\\V\\c*,D), (q ,q) G (dD x dD)\dG uniquely 
determines V and that k'(E,qo,q) given for all (qo,q) G (dD x dD)\dG uniquely 
determines V on D at fixed energy E > E'(\\V\\ C 2, D). 

2.3 Properties of the scattering operator. For equation (1.1) under condition (1.2b), 
the following is valid (see [Y]): for any (v-,x-) G B c x M, d , v- 7^ 0, the equation 
(1.1) under condition (1.2b) has a unique solution x G C 2 (M, M. d ) such that 

(2.4) x(t) = v-t + x- +y-(t), 

where y~(t) — > 0, y~(t) — >• 0, as t — > —00; in addition for almost any G 

S c x M d , t>_ ^ 0, 

(2.5) =v + t + x++y + (t), 

where v + 7^ 0, |i> + | < c, v+ = a(i>_,x_), x + = b(v-,x-), y+(t) — > 0, — > 
0, as t — > +00. 
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The map S : B c x R d ^ B c x R d given by the formulas 
(2.6) v+ = a(v-, X-), x+ = b(v-,X-) 

is called the scattering map for the equation (1.1) under condition (1.2b). The 
functions a(v-, X-), b(v-, x_) are called the scattering data for the equation (1.1) 
under condition (1.2b). 

By T>(S) we denote the domain of definition of S; by TZ(S) we denote the 
range of S (by definition, if (y-,x-) G T>(S), then v _ 7^ and a(v-,X-) 7^ 0). 

The map S has the following simple properties (see [Y]): for any (v,x) G 
B c x R d , (v, x) G V(S) if and only if (-v, x) G K(S); V(S) is an open set of B c x R d 
and Mes(( .B c x R d )\V(S)) = for the Lebesgue measure on B c x R d induced by 
the Lebesgue measure on R d x R d ; the map S : T>(S) — > 7Z(S) is continuous and 
preserves the element of volume; for any (v,x) G T>(S), a(v,x) 2 = v 2 . 

The map S restricted to 



E E = {(v-,x-) G B c x R d I |u_|= c yi- (jjj } 

is the scattering operator at fixed energy E and is denoted by Se- 

We will use the fact that the map S is uniquely determined by its restriction 
to M(S) = V(S) DM, where 

M = {(u_,x_) G B c x R d \v- ^ 0,v-x- = 0}. 

This observation is completely similar to the related observation of [N2], [Jl] and 
is based on the fact that if x(t) satisfies (1.1), then x(t + t ) also satisfies (1.1) for 
any to G R. In particular, the map S at fixed energy E is uniquely determined by 
its restriction to Me{S) = V(S) n Me, where Me = ^e riM. 

2.4 Inverse scattering problem and inverse boundary value problem. Assume that 

(2.7) VeC 2 (D,R). 

We consider equation (1.1) under condition (1.2a) and equation (1.1) under con- 
dition (1.2b). We shall connect the boundary value data k(E,q,qo) for E > 
E(\\V\\ C 2,D) and (q,qo) G (dD x dD)\dG, to the scattering data a, b. 

Proposition 2.1. Let E > E(\\V\\ C 2, D). Under condition (2.7), the follow- 
ing statement is valid: s(E, q , q), k(E, q , q) given for all (q, qo) G (dD x 8D)\dG, 
are determined uniquely by the scattering data a(v-,X-), b(v-,X-) given for all 
(v-,X-) G Me(S). The converse statement holds: s(E,qo,q), k(E,q ,q) given for 
all (q,qo) G (dD x 8D)\dG, determine uniquely the scattering data a(v-,X-), 
b(v-,X-) for all (v-,x-) G Me(S). 
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Proof of Proposition 2.1. First of all we introduce functions x, r_ and r + 
dependent on D. 

For (t>, x) G M n \{0} x IR n , x(^ 5 x) denotes the nonnegative number of points 
contained in the intersection of dD with the straight line parametrized by H. — > 
R n ,t I— > tt> + x. As .D is a strictly convex open subset of R n with C 2 boundary, 
X(u, x) < 2 for all t),iG IR n , i> 7^ 0. 

Let (v, x) G R n \{0} x R n . Assume that x(v, x) > 1. The real T-(v, x) denotes 
the smallest real number t such that T-(v,x)v + x G and the real t+(v,x) 
denotes the greatest real number t such that t+(v,x)v + x G dD (if x(v,x) = 1 
then T-(v, x) = t+(v, x)). 

Direct statement. Let (qo,q) G (dD x dD)\dG. Under conditions (2.7) and from 
(2.1) and (2.2), it follows that there exists a unique (v-,X-) G Me(S) such that 

= 2, 

g = + r_(v_, x_)v_, 

g = b(v-, X-) + t+(cl(v-, &(i>-, x_))a(f _, 

In addition, s(i£, g 0) <?) = T +( v -i x -) ~ r -( v -i x -) an d k(E, qo, q) = a(t>_, x_). 
Converse statement. Let (i>_,x_) G .M.e(S'). Under conditions (2.7), if x(t>_,x_) 
< 1 then (a(p-,X-),b(p-,X-)) = (p_,x_). 
Assume that x( v -i x -) = 2. Let 

q = x- + t-(v-,X-)v-. 

From (2.1) and (2.2) it follows that there is one and only one solution of the 
equation 

(2.8) -k(E, q, q ) = v-, q G 8D, q ^ q . 

We denote by q(v-,X-) the unique solution of (2.8). Hence we obtain 

a(v-,x-) =k(E,q ,q(y-,X-)), 

b(v-,x-) =q(v-,x-) - k(E,q ,q(v-,X-)) 

x (s(E,q ,q(v-,X-)) + r-(v-,x-)). 

Proposition 2.1 is proved. □ 

For a more complete discussion about connection between scattering data and 
boundary value data, see [N2] considering the non relativistic Newton equation. 



3. Inverse boundary value problem. 

In this Section, Problem 1 of Introduction is studied. Following [GN], we reduce 
the inverse boundary value problem to the problem of determining an isotropic 
Riemannian metric from its hodograph. 
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3.1 Hamiltonian system. Let E > c 2 + sup x6 £, V(x). Take V G C 2 (R n ,R) such 
that V\£) = V. We shall still denote V by V. Take an open neighborhood Q of D 
such that E 1 > c 2 + sup x6n V(x). The equation (1.1) in O is the Euler-Lagrange 

equation for the Lagrangian L defined by L(x,x) = —c 2 \Jl — — V(x), x G B c 
and a: G O. The Hamiltonian H associated to the Lagrangian L by Legendre's 

transform (with respect to x) is H(p,x) = c 2 \Jl + ^ + V(x) where p G lR n and 
a; G O. Then equation (1.1) in O is equivalent to the Hamilton's equation 

x = — (p,x), 

(3-D % H 

dH 

3.2 Maupertuis's principle. In this subsection we apply the Maupertuis's principle 
to the Hamiltonian system (3.1). 

Let (p(t),x(t)), t G [£1,^2], be a solution of (3.1). Let ^(t) = (p(t),x(t),t), 
t G [t u t 2 \. Then 

7 is a critical point of the functional J defined by 

= / pdx — H(p, x)dt on the set of the C 1 functions 

Jy 

7': [ti,t 2 H»xflx [t lt t 2 ], (p'(t),x'(t),t) 
with boundary conditions x'(ti) = x{t\) and x'fe) = #(£2)- 



(3.2) 



Let S denote the 2n — 1-dimensional smooth manifold {(p, x) G lR n x O | H(p, x) = 
E}. From (3.2), it follows that 

for any (p(t),x(t)), t G [£1,^2], solution of (3.1) with energy E 
and for any strictly increasing C 1 function <p from some closed interval 
[£-,£+] of R onto [£1,^2], the C 1 map 7 defined by 
(3 3) = W^W)' * e i s a critical point for the functional 

J defined by ^(7') = / prfx on the set of the C 1 functions 

Jy 

j' : [£_,£+] — > E, 1 1— > (p (t) , x' (t)) with boundary conditions 
x'(t-) = x{t\) and x'(£+) = x(t 2 )- 

Let y G C 2 ([ti, £2], ^) be such that y(t) 7^ 0, £ G [ti,^]- Let y be the strictly 
increasing C 1 function from [£i,£+] (£+ > 0) onto [£1,^2] defined by 4> y (ti) = £1 
and H(§:(<j)(t)y(<P(t)),y((P(t))),y((P(t))) = E, £ G [£i,£+], i.e. y is the function 



1 ^ E-V{ y (4> y (t))) ^_ 9 



t G [t!,t+], with initial datum ^(h) = h. Let j(t) = (f§ (0(0-* (*))£(*), y(t)), 
y(t)), t G [hM- Then, Jtf(t)) = -^f=|^)^(t))|rft. Hence, using 



which satisfies the ordinary differential equation 4> y (t) = c , ^ 

»t+ \<t>{t)y{4>{t)) 

that H(~f(t)) = E,t G [£1,^2], we obtain that 

j(7(0) = [ t \ v My(Hmmt))m)dt 
Jti 

(3.4) = f\ v , E {y{t))\y{t)\dt 



where ry iE (x) = cy ( g J (x) j — 1, x G O. 

From (3.3) and (3.4), it follows that if x(t), t G [ti , t 2 ] , is a solution of 
(1.1) in O with energy E, then x(t) is a critical point of the functional l{y) = 
ftf r v,E(y(t))\y(t)\dt defined on the set of the functions y G C 1 ([ti, t 2 ], ft) with 
boundary conditions y(ti) — x(ti) and yfo) — xfo) (Maupertuis's principle). As 
l(y) is the Riemannian length of the curve parametrized by y G C ([ti, £2], ^) for 
the Riemannian metric ry^ (a;) m ^ 5 one obtains that if x(t), t G [ti,^], is 
a solution of (1.1) with energy E, then x(£) composed with its parametrization 
by arclength (for the Riemannian metric rv,E(x)\dx\ in O) gives a geodesic of the 
Riemannian metric rv,E(x)\dx\ in O. 

For any solution x : [0, £+] — > Q of equation (1.1) in O with energy i?, the 
parametrization by arclength of x(t) is given by the strictly increasing C 2 function 
ip x from [0,t' + ] (t' + > 0) onto [0,£+] defined by the ordinary differential equation 

= /^-VWW.))^ With initial datum ^x(O) = °- 
I 3 ) _1 

Applying Maupertuis's principle we obtained the following Lemma. 

Lemma 3.1. Under the assumption E > c 2 + sup x6f2 V(x) the following 
statement is valid: for any solution x : [0, £+] — > O of equation (1.1) in O witt 
energy E, the map y : [0,t' + ] — > O defined by y(t) = x(ijj x (t)), t G [0, is a 
geodesic of the Riemannian metric rv,E(y)\dy\ in O which satisfies rv,E(y)\y\ = 1- 

We obtain, in particular, that trajectories {x(t)} of the multidimensional 
relativistic Newton equation in O with energy E coincide with the geodesies of 
Riemannian metric rv,E(x)\dx\ in O where \dx\ is the canonical euclidean metric 
on Q. (In connection with the Maupertuis's principle and analog of Lemma 3.1 
for the Newton equation in the nonrelativistic case, see for example Section 45 of 
[A].) 

3.3 Simple metrics. We recall the definition of a simple metric g in a bounded open 
subset U of M n with C 2 boundary (denoted by dU) (see for example [SU]). 
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Let U be a bounded open subset of R n with C 2 boundary (denoted by dU) 
and let g be a C 2 Riemannian metric in U. For x G dU the second fundamental 
form IT (with respect to g) of the boundary at x is defined on the tangent space 
T x (dU) of dU at x by the formula 

n(C)=s*(v c iv(£),o 

where C G T x (dU) and iV(x) denotes the unit outward normal vector to the bound- 
ary at x (g x (N(x), N(x)) = 1), and where ViV denotes the covariant derivative of 
the vector field N with respect to the Levi-Civita connection of the metric g. 

We say that g is simple in U, if the second fundamental form is positive 
definite at every point x G dU and every two points x, y G U are joint by an unique 
geodesic smoothly depending on x and y. The latter means that the mapping 
exp x : exp^iJJ) C T X U — > [7 is a diffeomorphism for any x E U, where exp x (v) 
denotes the point which is reached at time 1 by the geodesic in U which starts at 
x with the velocity v at time (T X U denotes the tangent space of U at the point 
x). 

As it was mentioned in [SU], if a Riemannian metric g is close enough to a 
fixed simple metric go in C 2 (U), then g is also simple. 

Here, as D is assumed to be a bounded strictly convex open subset of M. n 
with C 2 boundary, it follows that the euclidean metric \dx\ is simple in D. Hence, 
from the fact mentioned in [SU], it follows that there exists £7(||y||c2, D) such 

that for E > E(\\V\\ C *,D) the metric crv ^ (x) \dx\ = c 2 ^ ^P^ j - ^\dx\ is 

also simple in D. 

Hence for E > E(\\V\\c^, D) the metric r VyE (x)\dx\ is simple in D. 

Then one can consider properties (2.1) and (2.2) as consequences of Lemma 
3.1 and the fact that the metric rv,E(x)\dx\ is simple in D. 

Let ly,E denote the distance on D induced by the Riemannian metric 
^v,e(^)|^|- 

3.4 Properties of ly,E at fixed and sufficiently large energy E. Let E > E(\\V\\c2, 
D). From properties (2.1) and (2.2) (or from the fact that ry^E{x)\dx\ is simple), 
it follows that 

(3.5) ly : E G C(D x D,M), 

(3.6) ly E eC 2 ((DxD)\G,R), 

(3-7) maxd^^^U^^^D^d, 

/ q q \ | d 2 ly E s, . C 2 

( 3 - 8 ) I a „ (C>g)l < 



9Ci^j ' IC - x\' 

for (C,x) G (-D x -D)\G, C = (Cij ■■.(<;). s = (zi> ",a:d), and i = l..n, j = l..n, and 
where C\ and C2 depend on V and the map vy : E '■ dD x D -> S 71 " 1 , defined 
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by 



has the following properties: 

(3.10a) ^v,e £ C l (dD x _D,S n_1 ), 



the map : 9-D -»■ S n \ C -»• vv,e((, x) 

is a C 1 orientation preserving diffeomorphism from dD onto S r 



for x G D (where we choose the canonical orientation of S n 1 and the orientation of 
dD given by the canonical orientation of M n and the unit outward normal vector), 



k (E,x,() 



(3.10c) 



k (E 7 x,C)\ 
k(E,(,x) 



\k(E,(,x)y 

for (C,x) EdDx D. Note that from (2.3), (3.9) and (3.10c) one obtains 

(3.11) (^( C , x ),..,§i£ K , x )) W<-*> 



dx x ' ' dx n ' / 1 fc(g,c^Q 



for (C, x) edDx D. 

3.5 Determination of an isotropic Riemannian metric. We consider the following 
geometrical problem : 

at fixed energy E > E(\\V\\ C 2, D), does I v ,e(C,x), given for all ((,x) G dD x <9D, 
determine uniquely ry,E on D 7 

Muhometov- Romanov [MR], Beylkin [B] and Bernstein-Gerver [BG] study the 
question of determining an isotropic Riemannian metric from its hodograph. Re- 
sults in [B] and [BG] are obtained with smoothness conditions that are too strong 
so that one could apply these results to our problem. Therefore, for sake of con- 
sistency, we give results (Lemma 3.2 and Theorem 3.1) that already appear with 
stronger smoothness conditions in [B] and [BG]. 

We denote by uj ,v the n — 1 differential form on dD x D obtained in the 
following manner: 

- for x G D, let u)v,x be the pull-back of ujq by i/v,e,x where lvq denotes the 
canonical orientation form on § n_1 (i.e. uq(C,)(v i, ..,v n -±) = det(C, i>i, ••, v n -i), for 
C G § n_1 and v\, .., i> n _i G T^S n_1 ), 
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- for (C,x) G dD x D and for any (vi, .., v n -i) G T^ x ^(dD x D), 
u ,v{(, x)(v u .., v n _i) = wv,x(0(<r[c,x)M> ^(C^K-i)), 

where a : dD x D — > <9.D, l_ ^ C'> an d °"(^ x ) denotes the derivative (linear 
part) of a at (£, x). 

From smoothness of vv,e, o and u>o, it follows that a;o,v is a continuous n — 1 
form on dD x D. 

Now let A G M + and V U V 2 E C 2 (D,R) such that max(|| V x || C a , H^Hc 2 ) < A. 
Let E 1 > E(\,D). 

Consider the differential forms $ on (dD x dD)\G and $1 on (<9D x -D)\G 
defined by 

$o(C,z) = - (-^)~ s ~d x (l Vii E - 1 Vi ,e){(,x) A d c (lv 2:E - 1 Vi ,e){(,x) 

(3.12) A (dd c l Vl ,E(Cx)rA(dd c l V2 , E (Cx)) q i 

p+q=n—2 

for (C, x) G (d£> x <9L>)\<5, where d = d ( + d x , 

$!(C,x) = - (-1)^ [d^sCC,*) A (rfd C ;y 1 ,E((,x)) n - 1 

+ d x l V2! E{(, x) A (dd c lv 2 ,E((, ^)) n_1 - d x l Vl ,E((, x) A (dd c l V2 ,E(C, x)Y~ 

(3.13) -^Zv^bCCx) A (dd c / Vll B(C^)) n_1 ] , 

for (C, x) G (dD x D)\G, where d = d^ + d x . 

From (3.6), (3.7) and (3.8), it follows that $ is continuous on (dD x 
and integrable on dD x dD and $1 is continuous on (dD x D)\G and integrable 
on <9D x I). 

Lemma 3.2. Let Ag1 + an<iE> E(X,D). Let V U V 2 G C 2 (L>,K) such that 
max(||Vi||c 2 ) II Vz ||c 2 ) < A. The following equalities are valid: 

(3.14) / $ = / $1; 

JdDxdD JdDxD 
^-^yy$l(C,x) =(ry 1)E (x) n W ,V- 1 (C,^) +n/ 2 ,£(x) n W ,V2(C ;r ) 

(3.15) - V a .Zv 1)B (C,x)V a .Zv 2) B(C,a;) 

x (r VuE (x) n - 2 u 0!Vl ((,x) + r V2jE (x) n - 2 u 0! v 2 ((,x))) 
A (ixi A .. A e?x n , 

for (C, x) G x D, where V x l Vi , E ((, *) = (%r (C, x), .., %^(C, *)) for (£ x) G 
dD x D and i = 1,2. 

Equality (3.14) follows from regularization and Stokes' formula. Using Lemma 
3.2, we obtain the following Theorem of uniqueness and stability. 

12 



Theorem 3.1. Let A G R + and E > E(X,D). Let V U V 2 G C 2 (D,R) such 
that max(||Vi||cf2, ||V2||c 2 ) < A. The following estimate is valid: 

/ (r VuE (x) - r V2jE (x)) (rv^Eixy- 1 - ry^x)" -1 ) dx < 
Jd 

(3-16) J, (f) / $o. 

27T2(n-l)! J dDxdD 



Note that for Vi, V2 € C 3 and <9-D G C°° Theorem 3.1 follows directly from 
the stability estimate of [B] and [BG] for the problem of determining an isotropic 
Riemannian metric from its hodographs. Similar Remarks are also valid for Lemma 
3.2. For Vi,V2 G C 2 and <9.D G C 2 an estimate similar to (3.16) follows directly 
from a stability estimate of [MR] for the problem of determining an isotropic 
Riemannian metric from its hodographs. As we have followed Gerver-Nadirashvili's 
framework [GN] , we have chosen to extend the related results obtained namely in 
[BG] and [B] to the case of less smooth metrics. 
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